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High  tide; 
Low  tide; 
Flood  tide; 
Ebb  tide; 


Terminology 

the  instant  at  which  the  water  level  in  the  estuary 
attains  its  maximum  height; 

the  instant  at  which  the  water  level  in  the  estuary 
attains  its  minimum  height; 

the  time  interval  beginning  at  low  tide  and  ending 
at  high  tide; 

the  time  interval  beginning  at  high  tide  and  ending 
at  low  tide. 
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Introduction*  _ 

The  distribution  of  salt  or  any  other  solute  in  an 
estuary  Is  governed  by  the  dynamic  flow  pattern  and  the  extent 
to  which  portions  of  water  having  different  solute  concentrations 
nix  vith  each  other.  The  flow  and  the  mixing  are  the  result  of 
the  complex  interactions, of  many  factors,  e, g,  , the  flow  due  to 
the  river,  the  tidal  flow,  the  ground  seepage,  evaporation,  wind 
stresses,  gravity  flow,  etc.  In  many  estuaries  some  of  these 
factors  are  of  much  greater  importance  than  others.  In  the  follow- 
ing investigation  we  restrict  ourselves  to  estuaries  in  which  the 
tides  and  the  river  flow  are  the  predominant  agents,  and  we 
neglect  all  other  influences. 

We  furthermore  confine  our  attention  to  estuaries  which 
are  "essentially  vertically  homogeneous’1.  Ey  this  we  mean  that 
if  measurements  are  made  of  the  concentration  of  the  solute  by 
taking  samples  over  an  entire  transverse  cross-section  of  the 
estuary  (i.  e. , a vertical  ^section  which  is  perpendicular  "to  the 
direction  of  the  river  flow)  at  any  time  during  one  or  more  tidal 
cycles  (during  which  time  the  conditions  of  river  flow  and  tidal 
amplitude  are  assumed  to  be  unchanged),  then  the  values  of  the 
concentration  so  obtained  will  differ  from  one  another  by  a 
smaller  order  of  magnitude  than  will  the  values  obtained  by  taking 
measurements  at  various  positions  along  the  length  of  the  estuary, 
but  anywhere  in  a given  cross-section  and  at  any  time.  A typical 
case  of  an  "estuary  which  is  not  vertically  homogeneous  is  one • 
which  exhibits  two  substantially  distinct  layers  of  different 
salinities,  the  more  saline,  and  hence  denser  layer,  being  on 
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top.  In  such  an  estuary  the  differences  in  salinity  between  the 
two  layers  at  any  given  position  along  the  estuary  and  at  any 
time  will  be  greater  than  the  differences' between,  the  salinities 
at  points  in  the  same  layer  but  at  various  stations  along  the 
estuary.  The  division  of  estuaries  into  vertically  homogeneous 
and-  "stratified”  ones  is,  of  course,  not  a sharp  one,  and  cases 
which  do  not  clearly  fall  into  one  or  the  other  category  will  be 
encountered  in  practice. 

The  long  range  aim  of  a theoretical  study  of  estuaries 
is  to  permit  us  to  predict  the  distribution  of  any  solute  within 
the  estuary  x’rom  a knowledge  of  the s river  flow  and  the  tides, 
possibly  including  the  introduction  into  the  estuary,  or  the 
removal  therefrom,  of  any  solute  by  means  of  an  external  source 
or  sink,  such  as  the"  discharge  of  a pollutant  by  an  industrial 
plant.  A complete  theory  would  require  a study  of  the  detailed 
velocity  distribution  and  of  the  mixing -processes.  At  our  pres- 
ent stage  of  knowledge  ttie  development  of  such  a theory  is  much 
too  difficult.  We  must  content  ourselves  with  a more  macroscopic 
approach,  i.e.,  one  which  aims  only  at  predicting  some  convenient 
average  varlati on  of  solute  concentration  as  a function  of  posi- 
tion along  the  estuary  and  possibly  of  time.  Such  an  approach 
has  been  made  by  Ketchum  [I]*  and  by  Arons  and  Stommel  [2], 
Ketchum  subdivides  the  estuary  along  its  axis  into' a number  of 
volumes,  each  of  which  is  as  long  as  the  average  excursion  of  the 
tide  in  the  neighborhood  of  the  volume,  and  he  assumes  that  com- 
plete mixing  tabes  place  within,  each  volume  at  high  tide. 

* Numbers  in  square  brackets  refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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Ketehum's  method  is  essentially  a refinement  of  an  older  method 
which  assumes  that  all  the  water  in  the  estuary  mixes  completely 
at  high  tide*  Arons  and  Sterne 1 formulated  Ketchum's  idea  in 
terms  of  a diffusion  equation  which  they  assume  to  govern  the 
solute  distribution.  The  equation  involves  an  "eddy  diffusivity" 
whose  value  is  assumed  to  depend  on  the  amplitude  of  the  tidal 
velocity  and  the  tidal  excursion* 

The  following  investigation  represents  an  attempt  to 
derive  (rather  than  to  postulate}  diffusion  equations  by  assuming 
that  the ^essence  of  the  mixing  process  can  be  described  in  terms 
of  very  simple  physical  models*  It  is  hoped  that  this  approach 
will  give  new  insight  into  the  problem  by  showing  clearly  what 
tjhe  assumptions  are  that  lead  to  a certain  diffusion  equation, 

And  how  modifications  or  generalisations  of  the  assumptions  change 
the  equation;^  Moreover,  such  an  approach  may  provide  some  guid- 
ance in  the  development  of  diffusion  equations  -describing  other 
physical  models  in  which  various  mixing  processes  may  be  postu- 
lated* Solutions  of  the se  equations  may  then  be  dnmparedhrfith 
observation  or  experiment  to  determine  in  which  estuaries,  if  any, 
the  solute  distribution  is  adequately  described  by  one  or  another 
equation.  In  this  way  we  may  hope  to  moke  the  initial  steps  in 
formulating  laws  which  describe  the  solute  diffusion  in  various 
estuaries  and  in  beginning  to  understand  the  function  of  the 
various  physical  factors  involved, 

Exceptii^'  the  final  section  of.  this  report,  it  will  al- 
ways be  assumed  that  all  ..quantities  ai’e- periodic  in  time  with  the 
period  of.  a tidal  cycle.  In  particular,  this  moon-  f ho  , "the  time 
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average  of  the  solute  distribution  over  a tidal  cycle  is  the  same 
for  every  tidal  cycle. 

For  brevity  we  shall  usually  speak  o’  the  salinity,  but 
it  is  to  be  understood  that  the  investigation  applies  equally 
well  to  any  solute. 

We  begin  our  considerations  of  mixing  by  examining  an 
estuary  or  channel  in  which  the  motion  of  fluid  and  the  distribu- 
tion of  salinity  (or  the  distribution  of  any  pollutant  that  is 
introduced  into  the  estuary)  will  be  assumed  to  be  governed  by 
three  .factors: 

1)  The  tides,  which  are  assumed  here  to  cause  a periodic 
displacement  of  part  of  the  water  of  the  estuary  relative  to  the 
rest  of  the  water- and  hence  permit  the  possibility  of  the  mixing 

of  water  from  different  parts  of  the  estuary.  (Thus,  specifically, 
in  these  considerations  in  which  we  investigate  the  salinity  dis- 
tribution, this  tidal  motion  permits  a mixing  of  the  more  saline 
water  downstream  with  that  which  is  less  saline  further  upstream.) 
Our  basic  picture  will  assume  that  the  velocities  due  to  the 
river  are  much  smaller  than  those  due  to  the  tides  so  that  a 
pronounced  upstream  flow  occurs  during  flood  tide. 

2)  A small  scale  mixing  or  eddy  diffusion,  distances  moved 
by  the  fluid  because  of  this  mixing  mechanism  being  assumed  small 
relative  to  the  distances  which  the  fluid  moves  due  directly  to 
the  tides.  We  shall  not  inquire  into  the  detailed  dynamics  that 
cause  this  small  scale  mixing  but  will  concern  ourselves  only  with 
the  effect  of  the  mixing. 
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3)  The  river  which  produces  a net  flow  of  water  (and  hence 
also  of  salt)  downstream.  The  transport  of  water  due  to  the 
tides  and  subsequent  small  scale  nixing  will  be  found  to  produce 
a net  flux  of  salt  upstream  by  mixing  the  more  saline  water 
nearer  the  mouth  of  the  estuary  with  that  which  is  less  saline 
further  upstream.  These  two  net  fluxes  will  be  set  equal  in 
treating  the  equilibrium  condition  in  which  there  is  no  net  flux 
of  salt  upstream  over  a tidal  cycle. 

The  dynamics  which  are  due  to  the  tides  and  the  effect 
of  the  small  scale  mixing  as  well  as  to  the  river  are  hypothe- 
sised. The  salinity  distribution  then  follows  from  the  equations 
of  continuity  of  salt  and  fluid,  since  we  assume  that  the  net 
flux  of  salt  across  any  section  of  the  estuary  is  sero  during  a 
tidal  period.  In  the  derivation  of  the  salinity  dj  stribution  we 
shall  first  find  the  flux  that  would  be  produced  by  the  tides  and 
small  scale  mixing  in  the  absence  of  a river  flow#  When  next 
taking  the  river  flow  into  account  ve  assume  that  the  velocities 
which  it  causes  are  small  in  magnitude  when  compared  with  the 
velocities  caused  by  the  tides  during  most  of  the  tidal  cycle, 
and  hence  that  the  basic  picture  first  described  in  the  absence 


of  river  flow  may  be  preserved  when  the  river  flow  is  included. 


^ a u a n 

l IV?  V.  UC*1  li  x 


4 „ 4 ~ 4 


X"#  TA 

KsYI'J 


regions;  a lower  one  of  cross-sectional  area  b-^(x,t),  salinity 
Si ( x « t ) (dimensions  ML  in  which  the  water  moves  with  velocity 


u^(x,t)  and  an  upper  one  with  cross-sectional  area,  salinity  and 
velocity  b2(x,t),  S2(x,t),  ^(x,!)  respectively;  x being  the  dis- 
tance coordinate  along  the  length  of  the  channel  and  u-,  , U2  being 
the  velocities  due  to  the  tides. 
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the  channel  contains  fluid  of  two  kinds,  one  with  salinity  s^. 
moving  with  velocity  u-j_,  the  other  with  salinity  s2,  moving  with 
velocity  U2,  and  finally  that  the  total  cross-sectional  area  of 
the  fluid  of  salinity  s-j_  be  b^,  that  of  salinity  s2  be  b2*  For 
convenience,  however,  we  shall  speak  of  the  channel  as  if  it 
contained  two  distinct  layers* 


Kinetic  Considerations* 

We  begin  by  assuming  the  channel  to  be  of  constant 
width,  but  shall  generalize  later  to  include  the  case  in  which 
the  width  varies.  For  our  consideration  we  shall  need  the  equa- 
tions of  continuity  for  fluid  and  for  salt  and  some  assumptions 
about  b,  and  b2.  Further,  we  will  need  to  make  some  assumption 
regarding  the  transfer  of  salt  between  the  two  layers.  We  nave, 
then,  for  the  continuity  of  fl^id  in  the  lower  layer, 

3(b^Ui) 


5b, 

i 4. 

et  3x 


= 0 


(l) 


and  for  continuity  of  fluid  in  the  upper  layer 
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Although  we  assume  two  distinct  layers,  each  possessing  its  own 

velocity,  an  exchange  of  salt  between  the  two  layers  is  permitted 

so  that  the  equation  for  the  continuity  of  salt  is  for  the  total 

channel:  . . . 

dCbnS-.  + b0s0)  3 ( b-,  s i u,  + b0s0u0; 

i-o 2_2_2_  = (3) 


3t 


6x 


? 


Substituting  (1)  and  (2)  in  (3)  we  have 


-3St  esi,  u r^s2  9s2  i ~ 

[■— r=  + Ut  — — 1 + bo  f— ~ + u~  - — - 1 = 0. 

Lat  j-  ox  ^ -at  w ax 


(4) 


It  may  be  noted  that  in  writing  (1),  (2),  (3),  and  (4)  we  have 
assumed  that  the  fluid  in  each  layer  moves  with  uniform  velocity 
over  the  entire  depth  of  the  layer.  If  v/e  assume  instead  that 
the  fluid  motion  takes  place  over  depths  b-[,  bj,  (i.e.  , not  tho 
entire  depth  of  the  layers)  then  (1),  (2),  (3),  and  (4)  become, 
respectively, 


ai>, 

■*T 4- 

a(b-[  uT) 

= 0 

at 

ax 

*2  + 

a(b2  Up) 

= 0 

at 

dx 

a(b1  s, 

ww  J 

. + b2s2> 

^ 3(b{  slul  + b2  S2U2) 

at 

ax 

r u , t 9Si . , 3So  \ 0S„X 

,,  u-  — r*  + b u,  — -±)  + (bp  + b'  u-  £)  = 0. 

- at  1 1 ax  * at  2 ^ dx  ' 


(M) 


T n ,.rV»Q  4"  1 1 /\».i  m n4  n ■ 1 ;7-  --  x.  ■ - 1 » I 

««  n « wc  wxxx  qssuiuc,  ; lywe  vn  J.  , l/Ilcit  D, 


= D- 


. » 

Dp  = Dp. 


Before  considering  what  laws  might  be  formulated 


regarding  the  transfer  of  salt  between  the  two.  layers  we  shall 
consider  a spc  ^ f ic  simplified  model  which  will  help  in  the  formu- 
lation of  these  .avs  and  serve  as  a guide  in  the  solution  of  the 
problem  when  we  have  more  general  models.  We  assume,  for 


\J  hi  If 
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first  spec if io  case,  that 

U1  “ ® 
b,  = a 

1 (5) 

b2  = b + r\ 
r)  = C sin  ijt 


where  a,  b,  C and  w are  positive  constants.  The  equation 
fc2  = b + r]  assumes  that  the  tidal  wave-length  is  large  compared 
with  the  length  of  the  estuary  so  that  the  level  of  the  estuary 
rises  uniformly  over  its  length. 

Thus  we  have  a stationary  lower  layer  of  constant 
cross-sectional  area  and  an  upper  layer  the  area  of  which  varies 
-sinusoidally  with  frequency  where  2it/u)  = % is  the  tidal  period 

Assuming  thatrlrhe  fluid  in  the  upper  layer  moves  uniformly  over 
the  entire  area  b2  = b + rj,  we  may  use  the  continuity  equation 
(2)  to  determine 

9b2  9(b2u2)  + n)u2l 

at  ' ax  at  + ax 


a[(b  + q)u2] 

— Cu)  COS  wt  + — — — - ■■  ■'  — ■■  ■■  • 

dx 

Thus 

(b  + q )u2  = - £xu>  cos  wt 
or 

Kx  w cos  wt  / s \ 

u„  = - (6) 

b + C sin  wt 

where  we  have  set  the  arbitrary  function  of  t which  results' 
from  this  integration  equal  to  zero  sc  that  u^Ojt)  = 0 for  all  t 
We  keep  in  mind  that  the  region  of  validity  of  the  expressions 
for  q and  u2  probably  end  for  some  positive  value  of  x«  Equation 
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(6)  gives  the  velocity  u2(x,t)  of  a particle  in  terms  of  its 
position  x and  the  time  t.  The  position  of  the  particle,  l0tt) , 
as  function  of  time  may  be  found  by  substituting  1 2 £o:r  x and 
*J(t)  for  u2  in  (6)  and  integrating: 

w cos  ^t 


or 


so  that 


ch. 


2 _ 


dt 


b + C sin  wi 


d-t_  £u>  cos  u>t 

2 dt 

'C v-.  K -i_  J'  ^ * y-s  . » -f" 


b + C 


5-Lii  Ui  U 


log  ^2  = -*-°2  ^ (b  + C sin  ut) 


-1 


ln  = 


/ n \ 

C / / 


'2  (b  + £ sin  u)t) 
where  the  constant  of  integration  k may  be  determined  by  knowing 
the  position  at  one  particular  time: 


^2^o ^ ” 


b + £ sin  ut. 


from  which 


l2( t)  - 


l2(t0)(b  + C sin  wt0) 
b + X,  sin  wt 


(7’) 


Mixing  Considerations  — Calculation  oc’  the  Salt  Flux  * 

flaying  consider ea  the  motion  of  the  upper  layer 


caused  by  the  tides  directly.,  we  now  consider  the  second  factor 
which  governs  the  salinity  distribution,  the  small  scale  mixing 
which  produces  a transfer  of  salt  between  the  two  layers. 

V/e  assume  that  mixing  takes  place  Instantaneously 
at  high  and  low  tide  only,  and  further  that  the  mixing  at  each  of 


these  times  is  complete p i. e. , that  the  salinities  of  the  upper 


and  lower  layers  are  equal  following  each  of  these  mixings.  Me 
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shall  discuss  later  the  plausibility  of  this  assumption.  For  the 
present  ve  note  that  this  assumption  is  in  agreement  with  our 
concept  of  an  essentially  vertically  homogeneous  estaary. 

Our  object  is  to  calculate  the  flux  of  salt  past  a 
given  cross-section  during  a tidal  period  due  to  the  action  of 
the  tides.  Since  our  ultimate  goal  is  to  derive  a differential 
equation  for  the  salinity  as  a function  of  position  along  the 
estuary  we  shall  want  to  describe  the  flux  in  terms  of  the  salini- 
ty. The  differential  equation  will  express  the  fact  that  in  an 
estuary  in  which  conditions  do  not  change  from  one  tidal  cycle 
to  the  next  the  upstream  flux  of  salt  due  to  the  tides  must  be 
exactly  balanced  by  the  downstream  flux  due  to  the  river.  The 
flux  will  be  obtained  in  the  first  instance  (see  Eq.  (29))  from 
the  mixing  considerations  described  above  in  terms  of  the  salini- 
ties at  a fixed  position  but  at  different  instants  of  time, 
namely  those  just  ;jrior  to  and  immediately  after  the  mixing  at 
high  and  low  tide.  Using  the  equations  of  continuity  of  salt 
within  the  ur.per  layer  these  salinities  will  then  be  expressed 
in  te^ms  of  the  salinities  at  a fixed  time  and  at  various  x (see 


Eqs.  (30) , (33) , and  (34)).  From  these  a differential  equation 
for  the  salinity  as  a function  of  x is  derived. 

In  the  model  we  now  consider  the  salinity  of  the 
upper  layer,  as  measured  in  a frame  of  reference  moving  with  the 


velocity  of  the  upper  layer,  is  constant,  except  when  u>t  71/2, 
3r/2,  ...  and  hence  ..  . 


12  = 0 


r/ r , 3 75/  ^ , ...  . 


except  when  = 
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From  (8),  in  particular, 


(Vn  + l)it  - (4n  - Dm  , 

2(X5  — ~ s2^  + ^2 } 2uj u)> 


(9) 


n — 0 , 1 , + 2 j • • • 


where,  from  (7‘), 


x = SCbjJL) 

S2  b - C 


(10) 


That  is,  the  salinity  of  the  upper  layer  at  some 
point  x along  the  esbuary  just  before  mixing  atMg’n  tide  is 
equal  to  the  salinity  of  the  upper  layer  just  after  the  previous 
low  tide  mixing  at  the  point  x + further  downstream,  being 
the  distance  traveled  from  low  to  high  tide  by  that  section  of 
the  upper  layer  which  is  at  x at  high  tide.  In  similar  fashion. 


s20c,  lbu£l>i-  o)  =-  s2(x  - Kit  + O) 


where,  from  (71 ) 


2u 


(11) 


r»  — A -i_  T + P 

x*.  — '-'-J  a.  2 4.  j 


- filial* 


(12) 


b + c 

being  the  distance  traveled  from  high—' to  icw  tide  by  that  sec- 
tion of  the  upper  layer  which  is  at  x at  low  tide. 

Since  the  lower  layer  is  assumed  to  be  stationary, 

we  have  /).  _ . i \ . / ■ , _ %_ 

/ C-rn  + ih  i (H-n  - 

*-  — 35 0)  ■ sl<*>  — ^ 


0) 


(13) 


s,(x,  il-V-3—  - o)  = sx(x,  * 0),  (14) 


aw 


n = 0 , _f  1 , + 2 , 


• • ft  ft 


Further,  from  the  assumption  of  complete  mixing  after  high  and 
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low  tines,  we  have 


s,  (x, 


(4n  + 1 ) 7t 
2u> 


+ = So(x, 


(Ln  4-  i)n 


+ 0) 


CLUi 


!<*,  (>*  -3-^  4 0)  = Sa(x,  £^±-a>S  + 0) 


(lb') 

(16) 


n = 0,  + 1,  + 2,  . . . . 


At  this  point  we  leave  for  a moment  our  rather 
formalized  derivation  and  show  how  the  flux  over  a tidal  period 
may  be  derived  quite  simply  by  properly  neglecting  terms  of  rela- 
tive order  £/b.  We  shall  not  attempt  to  justify  the  particular 
approximations  that  are  made  therein  since  a more  rigorous  deriva- 
tion will  be  given  subsequently. 


Derivation  of  Diffusion  Equation  from  Elementary  Physical  Consider- 

.Slions. 

Let  us  begin  with  the  conditions  just  after  high 
tide  mixing  and  consider  the  transport  of  salt  during  a tidal 
cycle  across  the  section  at  x.  During  ebb  tide  a volume  of  water 
of  length  ^ ( see  Fig.  2)  and  cross-sectional  area  b + C b moves 
downstream  past  x.  Since  no  mixing  takes  place  the  volume  carries 
all  its  salt  content  past  the  section  x.  At  low  tide  it  instan- 
taneously mixes  with  the  more  saline  water  below  thus  acquiring 
a certain  amount  of  salt.  During  flood  tide  the  same  volume  of 
water  (now  of  cross-sectional  area  b - band  length  £ 2 ) 
moves  upstream  across  section  x.  At  high  tide  this  volume  occunies 
the  same  position  .as  it  did  at  the  beginning,  of  our  cycle,  but, 
having  acquired  an  amount  of  salt  at  lew  tide  It  has  a higher 
salinity  than  it  had  just  after  high  tide  mixing.  Moreover,  the 
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lower  layer  upstream  of  x,  having  lost  some  salt  during  low  tide 
mixing  with  a section  of  the  upper  layer  which  had  moved  down 
during  ebb  from  a position  further  upstream,  is  now  less  saline 
than  at  the  beginning  of  our  cycle.  Hence,  when  mixing  takes 
place  at  high  tide,  the  upper  layer  gives  up  salt  to  the  lower 
one.  The  cycle  is  now  complete  and  we  see  that  the  net  result 
has  been  'an  upstream  transport  or  flux  of  a certain  amount  of 
salt  past  the  section  x.  We  now  describe  this  process  in  mathe- 
matical terms. 


Let 


s2(x, 


(*+n  + 1)  it 

2w 


+ 0)  = Sjj(x) 


s2(x,  + 0)  = SL(x). 


(17) 

(18) 


For 


t «b,  K 


-V 


u2dt^ 


Cxw 


cos  ut  dt  = 


r> 


Jtc/2 


,t  m/2u 


from  (6),  Thus  we  may  call 


K = 


2xK 


(19) 


the  approximate  amplitude  of  tidal  displacement.  (Note  that  if  we 
neglect  terms  of  higher  order  in£  /b  in  (10)  and  (12),  we  have 
?1  * K •)  Then,  using  ( 11 ) , ( l4) , (16),  (17),  and  (18)  we 
may  write  the  equation  for  the  continuity  of  salt  at  low  tide  as 


aS^(x)  + bSH(x  - O 


= (a  -f  b)Sr(x) 


(20) 


and  that  at  high  tide  as 


aST(x)  + bST (x  + 0 = (a  4 b)3^(x) 


(21) 
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Fig.  2 

Schematic  diagram  showing  motion  of  water  due  to  tides 
and  mixing. 
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in  which  we  have  substituted  £ for  both  and  and  neglected 

terms  of  relative  order  £/b  (cf«  (20)  with  (33) )•  Now  the  salt 

transport  across  a section  due  to  the  tides  is  approximately 
^ , r 

b\  S^(x  + 'Od-t'w  b^Sjj(x  + p up,  during  flood  tide  and 


i 

bj  Sjj(x  - 'Od't''-'  b^Sp(x  - p down,  during  ebb  tide,  so  that  the 
net  flux  up  due  to  the  tides  during  a tidal  cycle  across  a section 


at  x is 


F = bUSL(x  + |)  w SgCx  - . 


(22) 


Returning  to  (20)  and  (21)  we  may,  for  C/b  « 1,  write  these  in 


the  form 


(X  + p + bSR(x  - |)  = (a  + b)SL(x  + |) 


( 20  * ) 


a^l/X  “ 2^  * bSI2x  + |)  - (a  ■+  b)SR(x  - |)  ( 21  * ) 

.respectively.  We  may  now  substitute  (201)  in  (22)  and  obtain  an 
equation  in  Sjj  alone: 


F - 


rb{sH(x  + r " bix  - |)3 


(23) 


* -rh  aH(x) 


plus  terms  of  higher  order  in  £,  and  hence  of  higher  order  in 


C/b.  This  result  is  identical  with  that  given  later  in  (3 7)  since 
2x  ^ 

from  (19)  £ = — — and  since  tlie  constant  c in  (37)  is  finally 
set  equal  to  re^o.  The  difference  in  sign  occurs  because  F is 


here  the  net  flux  upstream  whereas  in  (37)  F is  the  net  flux 


downstream. 


Similarly,  substituting  (21*)  in  (22)  we  may  obtain 
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an  equation  in  alone: 

ab£ 


v - 


[SL(X  + *)  - 


a + b 

s’(x) 

a + b L 


sl'x  -f» 


(21 


plus  terms  of  higher  order  in  Thus  if  we  call  s(x)  = 

SH(x)  + ST(x) 

- — the  average  salinity  during  a tidal  cycle,  then 

from  (23)  and  (24) 

F = S&IsL  (25) 

a + b 


Since  we  assume  that  there  is  no  net  increase  of 
salt  over  a tidal  cycle,  this  net  flux  upstream  must  eqvai  the 
net  flux  downstream  due  to  the  river  flow,  so  that 


r *s  = -aMi  is 
9 a + b dx 


(26) 


where  Rc  is  the  river  flow  per  tidal  cycle.  Equation  (26)  is  the 
integrated  form  of  the  diffusion  equation 


r0  4*  = 4-  ( 

u dx  dx 


ds) 


a + b dx 


(27) 


that  is  obtained  again  in  (41),  Using  (19), Eq,  (26)  may  be  inte 
grated,  as  is  done  in  Eqs.  (4l)  to  (43), 


Calculation  of  Salt  Flux,  continued . 

Although  this  direct  derivation  of  the  flux  and 
salinity  is  the  most  satisfactory  as  far  as  seeing  the  origin  of 
the  various  terms  in  (26)  is  concerned,  a more  formal  derivation 
in  which  the  terms  that  are  to.  be  neglected  are  clearly  exhibited 
was  found-  preferable  when  generalizing  the  various  factors  that 
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influence  the  final  salinity.  We  turn  therefore  to  the  formal 
computation  of  the  flux,  F(x),  of  salt  at  a fixed  section  x over 
a tidal  period  due  to  the  tides  and  small  scale  mixing  alone: 


2u)-0 

F(x)  = i bpUpS0  dt 

} C,  C.  w 

(28) 

W-n/2w  -0 

i 11 
11- 

(F  here  is  the  net  flux  downstream.)  In  the  integrand  (which  is 
a function  of  x and  t)  x is  held  constant  when  integrating  over 
t.  The  choice  of  the  limits  of  integration,  it/2w  - 0 to  57t/2w  - 0, 
is  arbitrary;  we  could  have  chosen  any  interval  of  length  2h:Aj. 
Since  we  want_to  speak  of  mixing  or  salt  transfer  between  the 
two  layers  "at"  high  or  low  tide,  i.e.  , when  ut  = --  or 

ut  = (n  = 0,  +1,  +2,  • • • ) it  is  convenient  to  have  these 

points  either  definitely  within  or  definitely  outside  of  the 
range  of  integration,  hence  the  "-0"  following  */2w  and  ?7i/2u  in 
the  limits  of  integration.  Since  the  integrand  in  (15)  is  finite 
for  all  t,  we  may  omit  the  intervals  rc/2u>  - 0 < t <~%/2u  + 0 
and  3,k/2w  - 0 < t < 3^/2w  + 0 and  write 

„(3n/2w-0  n5x/2u>-0 

t : 

F(x)  = b2u2s2dt  +|  b2u2s2dt.  (23 ') 

J V2u+0  J 3x/2u+0 
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Now 


gf  9(82^232) 
"dx  ax 

S ( 02^2  ^ 


- s2 


= - s 


3X 
9 b 

2 at 


b2u2 


5 S' 

w 


. 9S? 

4 g. 

2 ^ dx 


in  view  of  (2) 


_ ^ ras2  , „ ds2\  L 9^b2s2^ 

“ 62  (aT  u2  aF)  “it — • 


However,  in  either  of  the  intervals  ~-  + 0 <•  t ✓ -i2  - 0. 

— ~ 2u  1 


+ 0 <t  < - 0 we  have  -^2  - 


2u 


3x/ 2o^Q 


9t 


dx 


2u 

= 0,  so  that 


dF  - 
dx  " 


i 

J ti/2u+0 


a(b2s2 ) 


,5tc/2u-0 


at 


dt  - ! 


aCb-jSg) 


at 


dt 


J3^/2uh-0 


= (b  + Os2(x,  JL  + 0)  - Cb  - C)s2(x,  - 0)  — (29) 

+ (b  - C)s2(x,  + 0)  - (b  + C)s2(x,  I?  - 0) 

2u»  ^ 2w 


in  view  of  (5)* 


We  now  express  s2  at  t = — - 0,  + 0.  and  h.  ^oin 

zu  dxji  2ui 

terms  of  s?  at  t = -J-  + 0,  From  (11) 

S2(x,  ir  - 0)  = s2(x  - Kv  ■—  + 0)*  (30) 

In  order  to  express  s2  at  is  + 0 and  l£  - 0 in  te  ms  of  s0  at 
+ 0 we  need  the  equation  of  continuity  cf  salt  at  the  time  of 
mixing,  which  may  be  written  as 

as,(x.  - 0)  * (b  . C)S,b,  Aillih  . 0) 


= asitx 


2u 

(4n  + 1)tc 


,'U 


<?D 

0)  + (b  + C)s2(x,  llLtlh  + 0) 


-'3D 


Konr -5^20?/2 
at  high  tide,  and 


/„  (4n  + 3)u  , , ( (4n  + 3)71 

asn(Xj  g- ■ — — 0)  + (b  ~ ?3)sp(x,  ^ ~ 0) 

iU  £1U) 

= a31(x,  ttgJtUil  + 0)  + (b  - Os2(x,  iit^l-12*  + o) 


at  low  tide.  Substituting  (11),  (14),  and  (16)  in  (32)  we  have 


sa(x>  * ?>  = 

a i b~-g  tas2Cx,  f.  + 0)  + (b  - Os2(x  U-  + '0)j3|> 

and  from  (9)  and  (33)  we  have 

S2(x,  |s  - 0)  = 

a + b -1  [aS2U  + t’  it  + 0)  + (b  ' C)S2tX’  it  + 0,]C3lt) 

where-  it  is  to  be  noted  that  the  argument  of  the  term  multiplying 
(b  ) in  (34)  is  now  the  high  tide  position  of  that  section  of 
the  upper  layer  which  is  found  at  x + (2  at  low  tide,  i.e, , x,  and 
not  x + ^ Substituting  (30),  (33),  and  (34)  in  (29),  we 

have,  finally,  substituting  e = C/b 

s2^x  + 3533  + 

(35) 

S2(x  - r JL  + 0) ] 

£-UJ 

where,  from  (9')  and  from  (10*) 


dF  = [s^(Xj 

ax  a + b(l  - eT  l'  ’ 


ab(l  - e) 

! + b(i4  vy 


L So  ( X , 


oc  + K2 


x(l_  - e_) 
1 ’+  e 


♦ 


<35r) 
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Assuming  e « 1 we  expand  the  right  hand  side  of  (35’)  in  a povier 
series  in  e and  keep  only  the  terms  of  lowest  order  in  e,  which 
gives  „ 


dF  _ _ 


ab 


dx 


2 a + fc  dx 


2 ds2ix,  2w 


7“  + 0)  I 


dx 


(36) 


The  flux  of  salt  at  a section  at  x,  over  the  tidal  period 


is,  therefore. 


, ~2  . o 3s0(x.  J-  + 0) 

= . i&_  —a x2  _J — LS5 * c. 


2 a + b 


(37) 


it 

gi 

t-  a.  t u 0X 

D 

We  now  compute  the  flux  due  to  the  river  flow  (the 
third  factor  which  determines  the  salt  distribution)  and  add  this 
flux  to  the  flux  produced  by  the  tides  (given  in  (37)),  so  that 
t^ere  is  no  net  flux  of  salt  over  a tidal  period.  The  flux  due 
to  the  river  is  5V 2<*>-0 

G = Rs*(x, t)dt  (3.8| 

Jit/2u>-0 

where  R is  the  river  flow  across  any  section  and  has  the  dlmen- 
1 —1 

sions  of  LJT  , and  s*  is  related  to  s-,  and  s2  as  follows: 

If  the  river  produced  a flow  in  each  of  the  layers  pro- 
portional to  the  cross-section  of  the  layer,  then  we  might  define 
as-,  + (b  + t)  )s2 

s - — ± — ♦ If  the  river  acted  over  a cross-section 

a + b + r| 

Mx the  lower  B(x,t)  on  the 

'h;. . Bs p 

ucper  layer , then  we  would  have  s*  = ■ '-ry-  . However  , over 
a tidal  cycle  the  net  flux  produced  by  the  tides  and  sma l :l  scale 
mi.xing‘“ana  by  the  river  flow. is  r,  so  that  the  salinity.. of 
each  of  the  layers  is  periodic*  Thus  we  may  write,  for  any  time 

5x 


t ) < Sp!,x,  - G) , which,  together  yi  th  "f-j-CT '■  .133) , .( 3V) , 
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(9*)  and  (12)  give,  for  any  time  t, 

<tr  r 

s-^Xjt)  = s2(x,  + o)  + o(-) 

s2(x,t)  = s2(x,  JL  + 0)  + 0(— ) 
s*(xrt)  = — 


As-j_  + Bsp 


B ~ “ s2(x>  S + 0)  + 0(l)- 


(39) 


Thus  (38)  becomes 


G “ “ R [s2(x,  JL  + 0)  + 0(|)  ] . 


(1) 


(40) 


Derivation  of  Salinity. 

In  adding  the  flux  due  to  tides  and  small  scale  mixing 
to  that- due  to  the  river  flow  we  neglect  the  terms  in  (40)  of 
0(~),  Having  shown  in  (39)  that  the  time  variation  of  the  salin- 

b 1 

±ty,±n  either  layer  is  of  0(i) , we  may  again  neglect  terms  of 

¥ D 

O(-g)  and  replace  32(x,  JL  + 0)  In  (37)  arid  (40)  by  s(x),  which 
may  be  interpreted  to,  be  the  salinity  at  a to  within  terms  of 
relative  order  | at  any  time.  Adding  F and  G in  (37)  and  (40) 

and  setting  their  sum  equal  to  zero  we  then  have 

,2 


ur 


ab  2 

^2  a + b dx 


C , % R„ 


(41) 


the  solution  of  which  is 


u> 


- , . 2jtR 

2tcI  1oS  (“  + c;  = 


b2  a 


+ c 


ab 


+ b 


so  that 


O r*» 

2 IT  n 


u! 


s(x)  - - .4. o 

2-kH  2tiR 


Cf  — it 
0)  x 


(4l 1 ) 


where 
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k = 71  a b b r 

2 cj  ab  ^2 


<4l») 


We  now  set  c = 0 so  that  s(C)  = 0,  which  is  consistent  with  the 
fact  that  112(0, t)  = for  all  t so  that  no  salt  can  get  upstream 
of  the  section  at  x = 0. 

The  constant  c * is  determined  by  knowing  the  salinity 
at  some  particular  point: 


s(L)  = 


2iR 


e 


0) 


c “ 


so  that 


and  hence 


isS  c,  | 

e w = s(L.)e 


a — c 
2uR 


- 1) 

s(x)  = s(L)e  L x * 


(42) 


In  particular,  if  we  denote  the  ocean  salinity  by  d,  and  define 
the”length  L of  "the  estuary  to  be  the  distance  between  the  point 
at  which  s(x)  = o and  s(x)  = d,  then 


sU).  = 

d 


I a - -l) 


03) 


It  should  be  noted  that  Arons  and  Stommel  [2]  arrive  at  this  same 
result  using  dimensional  analysis  and  call  k/L  the  flushing  number. 
In  terms  of  the  notation  used  here  their  flushing  number  equals 

— v — where  B is  a constant  which  remains  undetermined  in  the 
2BtT  UL 

paper  of  Arons  and  Stommel.  It  arises  in  their  paper  from  the 
assumption  that  the  eddy  diffusivity  A of  the  salt  transfer  equa- 

’ Q 

tion  *at  + u "ax  = "ex  which  they  postulate  as  a basis  for 

their  considerations,  is  given  by  A = 2B£0U0  where  2£0  is  the  total 
excursion,  over  a tidal  period,  of  a particle  due  to  the  tides, 

UQ  is  the  amplitude  of  the  tidal  velocity  and  u is  the  river 
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velocity.  Equating  their  flushing  number  to  k/L  (see  Eq.  (k-l,!)), 
we  have 

B = JL  (¥f) 

Tib 

which  relates  B to  the  parameters  of  this  paper* 

Calculation  of  Flux  for  Generalized  Mixing  Process. 

Let  us  now  consider  the  various  ways  in  which  we  might 
make  the  present  model  less  restrictive,  in  order  that  it  have  a 
better  chance  of  corresponding  with  reality,  and  investigate  how 
the  salinity  distribution  is  affected  by  the  modifications  which 
are  introduced. 

We  first  consider  a model  in  which  the  dynamics  are 

those  previously  described  but  instead  of  assuming  mixing  to  take 

place  at  high  and  low  tide  alone  (t  = JL  and  t = 4^,  respectively,. 

2 u>  2 ur 

we  assume  that  there  is  mixing  at  t = JL,  JdL,  i*2  , Moreover, 

2u  <ji  2u>7  u>  1 

we  shall  not  demand  that  the  mixing  be  complete  at  each  of  these 
times,  but  rather  only  that  at  each  mixing  a certain  quantity  of 
fluid  from  the  top  layer  exchanges  places  with  an  equal  quantity 
of  fluid  from  the  bottom  layer.  Let  Q(x,t)  be  the  cross-sectional 
area  of  the  section  of  fluid  that  is  exchanged  during  mixing  at 
time  t,  and  position  x along  the  estuary.  For  brevity  let 


s ~ = s1(x,t  - 0)  s^+  = S^CXjt  4-  o) 

s2-  = S2(x,t  - 0)  s2+  = Sp(x,t  4-  C). 

Then  in  view  of  the  continuity  of  salt 


(M-5) 


( b-^  — W ) S-j  4-  Q s2  — b^s-^ 


(*+6) 


(T„  cf.  Op''  /O 
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whlch  may  be  written  as 


and 


- b1(s1+  - 3^“)  = Q(s}“  - s2") 
(b2  - Q)s2“  + Qs^“  = b2s2+ 


(46 ‘) 
(47) 


which  may  be  written  as 


b^(s2+  - s2“)  = Q( s2“ 


” s ^ ) • 


(47* ) 


If  the  exchange  o"  salt  between  the  two  layers  occurs  continuously 
rather  than  at  discrete  points  in  time,  then  (46*)  and  (47’)  read, 
respectively, 


- bi( 


QS^(x,t)  9s^(x,t) 


at 

aso(x,t)  aso(x.t) 

- b2( + u2  — ^ ) = Q[s2(x,t)  - s1(x,t)]  (47") 


+ u 


3x 

>(x 

"ax 


) = Qfs^Xjt)  - s2(x,t)]  (46") 


1« e« , in  either  case,  the  rate  of  discharge  of  salt  from  any  one 
layer  into  the  other  is  proportional  to  the  difference  of  the 
salinities  of  the  two  layers#  It  may  be  noted  that  (46")  and 
(47")  satisfy  Eq.  (4)  for  the  continuity  of  salt: 


>3®“}  . 3 S ^ 

b (”7*  + U.  — ) 4-  bo(r  v"*  4 U„ ^ ) = 0* 

l i ax  2 at  2 ax 


Reexamining  the  case  in  which  we  had  complete  mixing,  it  may  be 
noted,  substituting  Cl 5)  in  (31)  and  (16)  in  (32),  that  both  of 
these  equations  may  be  expressed  in  our  present  notation  by  the 

13  -4.^  p 

eauation  4 b2s2“  - (b-^  * b2)s2+  from  which  s ^ .^.3. 

so  that 


— b-.( 


J2 


- Sr 


1 A 


b1s1“  4 b0s2~l. 

- 1 -4^—  3 

-•  l s2~  - si'  .3 


(M3) 


tr 
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or,  alternately,  since  both  (31)  and  (32)  may  be  written  as 


Vl 


o_  >3  c 
’ °2S2 


b^Cs^  - sL")  + b2(s2+  - s2“)  - 0, 


(49*) 


we  may  also  write  (43)  in  the  forms 


- t>2.(  sl+  " sl“^  = 


TTX  Csi'  ' s2‘] 


(48*) 


Thus  for  complete  mixing  Q '=  — 1— r-  , from  which  we  have  both 

JL  + _L. 

bl  b2 

Q < b^  and  Q < b2,  as  might  be  hoped. 

In  postulating  a mixing  that  is  lessMdian  complete,  we 


shall  assume 


q(t) 

1 . JL  1 


(50) 


where  a will  be  permitted  to  oe  a periodic  function  of  t (with 
period  2%/u)  independent  of  x,  and  0 < a £ 1.  It  is  zero  at 
times  other  than  nit/2u  , i.e. , when  there  is  no  mixing. 

In  the  case  in  which  we  assumed  mixing  at  high  and  low 
tides  only  the  flux  was  given  by  (28)  and  (28*).  To  obtain  the 
equation  for  the  present  model  which  will  be  analogous  to  (28*) 
we  omit  the  intervals  about  the  four  points  nft/2w,  n = 1,2, 3, 4. 
Following  the  work  previously  done,  we  arrive  at  the  equation 
which  corresponds  to  (29)  and"  reads 

^ = (b  + C)s2^x>  'gw  + “ bs2(x,  J - 0) 

+ bS2(x>  Z+  0)  " (b  “ C)S2(X*  k * 0)  (51) 


+ (b  - C)s2(x,  + 0)  - bs2(x,  ~ - 0) 

+ bs5(x,  — + 0)  - (b  + ^)sp(x,  111  - 0). 

w W ci  1 j 
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As  before,  we  must  express  the  salinities  at  mt/2uj  + 0 
(n  = 2,3,4)  and  at  nx/2u>  - 0(n  = 2, 3,4, 5)  in  terms  of  the  salin- 
ity at  -it/2u  + 0;  this  time  through  the  use  of  (7*),  (46*),  (47?), 
(49)  and  (50).  We  now  proceed  to  do  this.  We  have  (for  the 
reasons  mentioned  in  connection  with  the  derivation  of  (9)  and 
(ID), 

s2(x,  § - 0)  = s2(x  - A1}  + 0)  (52) 


where , from  ( 7 ’ ) 


xb 

b + C * 


being  the  distance  travelled  from  high  to  mid-tide  by  that 
section  of  the  upper  layer  which  is  at  x at  mid-tide.  Next, 
from  (4  7' ) and  (50)  we  have 


-b[s2(x,  Zj  + 0)-s2(x,  A - 0)]  = _~2l_[s2(x,  g - 0)-s1(x,  21  - 091(54) 

“ T* 


~ ^ ’ 7 u .. 

* 7 u 

• j 

1 « 

r JL  “ w 

a b 

lihere 

a„ 

n 

= a ( t 

nrx  N 
~ 2u 

(54‘) 

Noting  that 

si(x> J - 

0)  = 

Sl(x,  £ + 0) 

(55) 

for  the  reasons  mentioned  in  connection  with  (13)  and  (14)  we 
have,  upon  substitution  of  (52)  and  (55)  in  (54), 


s?<x»  5 + °)  ~ 


* t1  - ^i]s2u  * Ai-  * + 0>+  it  si(x’  * * 0)- 


(56) 


Further,  similarly  to  (52),  we  have 


3U 


s2(x,  ~ - 0)  = s2(x  - A2,  L + 0) 


U) 


where,  from  (7‘), 


x - 


h. 


x(b  - C) 


(57) 

(58) 
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Ths  physical  meaning  of  Ap  as  well  as  all  further  A^  (i  =1, 

...  } 8)  which  will  be  defined  below  can  be  determined  from  Eq» 
(7f ). 

Substituting  (57)  in  (56)  we  have 


s2U>  “ E1 

where,  from  C 7’) 


act* 

g 

a+b 


4s 


2(x-A3, 


2w 


+0)  + 


^2 

a-<b 


It 

2u 


+ 0 ) 


b + C 


(59) 

(60) 


Next,  from  (4?*),  (50)  and  (54.»), 


- (b  - OCs2(x,  ||  + 0)  - s2(x,  - 0)]  = 


[s2(x,  £ - 0)  - Sl(x,  - 0)], 

a b-C 


(61) 


Now,  similarly  to  (55)  we  have 


3iCx>  IS  - °>  = si<*>  I * °) 

and  from  (46  1 ) end  ( 50) 

- a[sx(x,  2 + 0)  -s  (x,  & - 0)]  = T^r  [ s-i  (x.  21 

w 1 «*>  lil  1 'lil 

a + b 

sp(x,  ~ 

d.  7 U) 

v.iich  may  be  written  as 


(62) 


0)  - 
0)] 


(63) 


'lCx>  - * 0)  = [i  - ^]Sl(x,  | - c)  - s2(x,  a - 0) 


'=  t1  - iSiKCx,  -fe  + 0) 


ha- 


ul 


(63*) 


a+b^'”’  2^  T w T £3b  S2(x  - Al»  ^ + P) 
from  (55)  and  (52).  Thus  from  (63')  and  (62)  we  have 
3 Tt  „„  r - ^5:  , „■  ba. 


S1(X=  55- c)  = n -7  4 


^)Sl(x,  Ji  - 0)  + — ^ s?(x  - A , 4-  + 0) 
^ 2u  a+b  ^ i 2u 

(64) 
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and  substituting  (690  and  (59)  in  (6l)  we  have 


•a<*.  I * 0)  = ri  - Jgj  1 !C1  - s2(x  - V £ + 0)  ( 

J flrr_  ! 


aa 


! + a^S  sl(x  " A2»  * 

•j 

ba0  _ 

ru  - si(x>  & + 0) 


Next,  in  analogy  to  (52),  we  have 


/^TT 

s2(x,  — - - 0)  = s2(x  + A 


4’IS  + 0) 


where,  from  (7f), 


a xb 
x \ = FT-r  * 


Substituting  (66)  in  (65)  we  have,  using  (7*)., 


S2(*,  - - o)  - -C1  - -2^]  ^ 


‘ FrS1  32<x  * t5t  + 0) 


1>  2u 


31<X>  oS  + °> 


aJ-b  1 - 3 2W 


;■[!  - 3X(X  + A„  J2  * 0)  ■] 


a+b-C  <v 


U'-*  o , TT  . f 

+ s2(x+ V^  + 0>  i 


where 


A xo 

x 5 = -5 5 • 

5 ^ - C 


Now,  from  (9-7*),  (50)  and  (5^0, 


- b[s2(x,  ~ + 0)  - ag(x,.,|j  - 0)]  = 

a9-  f „ i ^ 2k  . _ / „ 2it  on 


1 1 

«*dfc  x *s» 

a b 


[ s2(x,  ~ - 0)  - s1Cx,  - 0)] 
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and  analogously  to  (55), 


Sl(x,  ^ - 0)  = Sl(x,  ^ + 0). 

Further,  from  (46’)  and  (50), 

- a tsr(x>  + °>  - si<x>  |S  - o)]  = 


^-T  tsi(x>  Is  - °)  - s2^x.  Is  - °n 


a b-C 


Ironi  which 


Bit 

S1CX’  2S  + 0)  = 


a . ihH!n  Sl(X,  g - 

a + b-£ 


..(b-C)a^  . . 

0)  + — --  — . B^SgCx, 


(71) 


(72) 


a + 


b-C 


- 0).  (73) 


Substituting  (59)  and  (64)  in  (73)  ve  have 

ba 

(b-C) a 


(x,  + 0)  = (X 


a-»-b-£ 


-h 


[!  ._Ig]  sl(x,  £+S» 


< 

I 


bcu 

+ t£  s2(x  - ai>  rs + 0)i 


| (1  - - — il  Sp(x  ■-  -S  + 0)j 
(b-C) a.  f a+b  1 2 3’  ! 


(74) 


a+b-C 


a a. 


^ + iTb  2l(x  “ ^2 9 + °^J  * 

Thus,  substituting  (74)  in  the  right  hand  side  of  (71)  and  the 
result  of  that  substitution  as  well  as  (68)  in  (70),  we  have 
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aa  r__.  aa^ 

’2'"*'  — + W = [1  - -^rg]  1 1-1 


r 


M*t  rr 


(75) 


! 

t 

b 

r\ 


a+b-£J  ■%  ado 


(1  - 
+ 


aa2 

a+b 


ls2(I'^!s; 40)) 


S.(x, 


71 


— — o ■=—  + 0)  }f 

a+b  1 7 2u 


aaa_  J 


it1  - ^t]31<*+\>2"  ^>|f 

ba2  u 


+ ZZbk  1 + ITb  S2U^5^  +0)j  I 


(b-C)a^ 

iCi  - -}i 


a+b  j 

*Y 


*■«  j V_ 

QTU* 


(b-t)a. 


C 1 + Xtl  s2!W,1>^5 

L 

" a CL 

Cl  - — gjs^Cx-^,^  40) 


1 


<1 

J 


L 


a+b-C  + 


+ S sl(x-A2>jE  40  J 


Finally,,  we  have  in  analogy  to  (52} 

i:  0)  '*  ^2U  :+  V%+ 

where,  from  (7*) 


x(b  + O 
X * %,:■=  -;■■■•- 


bstituting  <76)  in  (75)  we  have,  using  (71) 


% - o)  - [i 

c w 


aa 

+ 

a+b 


1 • 

an. 


• au.  1 


(78  r 


'<f 


aa^ 

i+b-r  ! 

^ -V 


+ 0) 

(77) 

?»} 

/ 

aa2  IS  (- 
a+b  jS2U?  2i 

+ 0) 

? 

I 

aau 

1 

1 

‘^+b  s1(x+A6» 

+ 

O 

W 

l 

ba0 

— fe  js,  (x+A  , 
a+b  1 7 

X-  - o)1 

2u>  l 

iy. 

boa, 

r.t0)  i 

—A  s (x+A,  , 

a+b  ^ + 

4 — *>-  T "» -i « " . 

j.  vi*  idEci 
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(H)«,  4 * ^f3si(x+i6>  ■& 


(?8)  COD.t* 


[1  - 


a+b-  £ 


U O *TT 

+ S2(X’  S + 0) 


(b-C)a  ;[l  ' i7b]s2(r~A2»  f'  +0)]f 


a+fc-C 


act  _ ( 

* Sl(x-A8>f.  +0)' 


where 


* *7 


x(b  4-  £) 

b - C 


X - * - 


2 2 

x(b  - C ) 


(79) 


We  now  substitute  (52),  (56),  (59),  (65),  (68),  (75),  and  (78) 

In  (51)  and,  making  use  of  (53),  (58),  (60),  (67),  (69),  (77), 
and  (79),  expand  the  resulting  expression  for  dF/dx  in  a power 
series  inC  /b,  keeping  only  the  terms  of  lowest  order  in  C/b. 

From  this  procedure  we  obtain  an  expression  involving  both 
s^(x,  n/2u)  + 0)  and  S2(x,  "n:/2w  + 0).  In  order  to  simplify  this 
expression  we  shall  make  the  additional  assumption  that  the  mixing 
at  high  tide  is  complete  (i.e.  , that  3-^  = 1 and  s2(x,  it/2u>  + 0)  = 
s^(x,  -rt/2u)  + 0)),  We  then  obtain 


j-n  y*2  p 5So(x,  *j~ ■ 4-  0) 

dr  _ _ 4 abM  5 r 2 * 2co  -1 

dx  , 2 a+b  fix  fix 


(80) 


M = + 4a^  + a-,  - 2a^a_  - a„3,  - 2a.  a.  + a a a, 

3 m-  2 j 24  34  234 


(81) 


It  may  further  be  noted  that  if  we  assumed  the  mixing  at  low  tide 
to  be  complete  ( =1)  and  that  at  high  tide  arbitrary  then  the 
above  equation  is  changed  only  in  that  must  be  replaced  by 
throughout.  Comparing  (80)  with  (36)  we  may  note  that  these  equa- 
tions are  identical  except  that  the  coefficient  !,4"  in  (36)  is 
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replaced  by  M which  reduces  to  4 If  we  let  cu,  = = 0 and  = 1, 

as  it  should.  It  may  also  be  noted  that  the  above  coefficient 
reduces  to  unity  if  we  set  0^=0^=  1,  = 0,  This  is  as  it 

should  be,  since  in  general  the  flux  is  proportional  to  the  square 
of  the  total  excursion  of.  the  upper  layer  during  a tidal  cycle,, 
and,  as  far  as  the  mixing  process  is  concerned,  the  upper 
layer  has  effectively  moved  (if  a2  = = 1,  = 0)  only  half 

the  distance  moved  in  the  first  model.  Following  the  work  pre- 
viously done  in  going  from  Eqs,  (.36)  to  (42),  we  may  note  that 
with  this  more  general  model  we  arrive  at  Eqs,  (42)  and  (43)  for 
the  salinity,,  except  that  the  parameter  k in  each  of  thesb  equa- 
tions must  be  replaced  by 


k*  = 


2*<a  + b)  jjag 


Id 


ab 


2 

XT  m 


In  (82}  it  should  be  noted  that  since  0 ^ ct^  = 2,3,4  we 

may  write 


.M  s 2a^>.  1 - d'g)  +-2CI2CI  - c^)  + + a^(l 


V 


i 0 

since  each  of  the  five 


(63) 


terms  in  the  second  expression  is  non- 


negative.  Noting ithat  we  may  also  write 

M = 4 - a2  - - 2(1  - a2>(l  - a^)  - 2(1  - a^Hl  - a^-d-a^a^ 


(84) 

since  each  of  the  five  terms  following  the  4 in  the  second  expres- 
sion is  negative,  we  have  from  (83)  ~*id  (84) 


C £ M ^ 4 
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where,  from  (83}  the  left  hand  equality  is-  -satisfied  only  if 

'=  0^  0 “and  the  right  hand  equality  is  satisfied  only 

if  ci2  = = 0 and  — 1*  If  a and  b can  be  obtained  from  obser- 

vations then  by  matching  the  observed  salinity  with  that  predicted 
by  this  theoretical  model,  it  should  be  possible  tcx  obtain  an 
estimate  of  M and  hence  to  obtain  an  idea  of  where  during  the 
tidal  cycle  most  of  the  'mixing  sakes  place,  large  values  of  M 
correspondingly  a state  in  which  most  of  the  mixing  takes  place 
at  high  and  at  low  tide. 


Extension  of  Analysis  to  Channel  of  Varying  Cross-Section. 

We  now  consider  a different  extension  of  our  original 
model.  As  first  assumed,  we  again  postulate  mixing  to  occur  only 
at  high  and  low  tides,  and  that  the  mixing  at  each  of  these  times 
is  complete.  The  dynamics  are  essentially  those  originally  des- 
cribed, However,  the  cross-sectional  area  of  "the  stationary  layer 
is  now  a(x)  in  place  of  the  constant  a,  the  average  area  of  the 
moving  layer  is  b(x),  and  the  time  varying  part  of  the  moving 
layer  is  £(x)  sin  wt  in  place  of  C sin  uft.  Following  the  calcu- 
lations connected  with  the  first  model,  we  see  that  Eq.  <29)  is 
changed  only  in  that  now  b = b(x),  C ~ C(x),  Equation  (30): 

3U 

sp(x,  ■—  - 0)  = Sp(x  - E , ~ + 0)  must  be  modified  only  in  that 
Z-l  is  no  longer  given  by  (12)  since  Eq.  (6),  from  v/hich  it  was 
derived,  must  be  modified  as  is  shown  in  (88).  With  this  modifi- 
cation implied,  (33)  must  now  read 


So 


(x, 


ill 

2u 


+ 0)  = 


a 


«)+b(x)-CR7  [a(x)s2<x>  ir 


2y 


0} 


f V f v \ 

\u\a| 


- C(x)  )s/x 


S-i  j 
i. 


2w 


* o)3 
(85) 
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Equation  (3^)  must  be  slightly  modified,  however,  since  substitut 
ing  (9)  in  (85)  (rather  than  in  (33))  we  have 

s2(*»  Is  - 0>  = 

[a(x+^2)s2(x+^2,  +0)+(b(x+C2)-  ((x^2))s2(x,  jL  + 0)3 

a(x+^2)  ■+  b(x+^2)-  C(x+^2) 


(86) 


+ W j 

Vi  -A 


in  which  it  must  be  recalled  that  $2  is  no  longer  given  by  (10)* 
Substituting  (30),  (85)  and  (86)  in  (29)  we  now  have 

dF  a(x)  [b(x)  + C(x)  ] r , n ^ „ * 

dx  a(x)+b(x)  - £(x)  ^ ^ ^ ^ 2u: 

a(x+C  ) [b(x)+C(x)  ] _ 

+ [Sp(x,  A.  +0)-s?(x+^„,  ~ + 0)3* 

aCx+^Wx+^-ttx+^r  2 2“  2 2’  2“ 

We  must  now  express  ^(x)  and  ^(x)  in  terms  of*  b(x)  and  C(x), 
and  hence  must  recalculate  Eq.  (7)  giving  the  position  of  a parti 
cle  as  a function  of  time  for  the  case  in  which  a,  b,  £ are 
variable* 

Returning  to  Eq*  (2)  for  the  continuity  of-fluid  in  ohe 
upper  layer  we  have 

9 [ ( b(x)  + £ (x)  sin  ut)u2]  - - -r|[b(x)  + C(x)sin  ut] 


dx 


- - £(x)u  cos  wt 


which,  upon  integrating,  gives 


[b(x)  + C(x)sin  wtlu2  = 


x 

n 


- u cos  ut  ! C(x’  )dx* 


V 


0 


♦See  remark  following  Eq.  (6), 


Ilonr-  56202/2 


35 


or 


JS. 

r 

u)  cos  uft , j 
0 


C(x* )dx‘ 


u = , 

b<x)  + C(x)  sin  wt 


(88) 


Following  the  procedure  used  in  the  derivation  of  (?’)  we  write 

7 

w cos  ut  C(x  ’ )dxf 
d£  _ JO 

b())  + r (i.)  sin  tot 


which  may.  be  written  in  the  form 

l 


[sin  <ot  j tCxOdx*]  = - b(£) 


(89) 


0 


regarding  t as  a function  of t.  From  (89),  recalling  that  x + ^ 
is  the  low  tide  position  of  a particle  which  is  at  x at  high  tide 
we  have 


l 


[sin  wt 


l't=x+^2?t=3x/2u 


0 


or 


C(x')dxf  i ^ 

J-teuc,  t=m/2w 

:^2 

6 


S+<2 


X 

x+K, 


b(l)dl 


r 


((x')dx*  +j  ((x')dx’  = j bCZ/dl, 
0 x 

Similarly,  we  have 


(90) 


t x 

{ 1 ^x,t=37t/2uj  f 

[sin  tot  C(x')dxf  = - \ M-Od-i 

\j  %PX“  ^ , t=  m/2  u>  vj 


or 


0 


0 


C(x’)dx’  + 


: r/2 t 

x 

f' 


X-£, 


((x  * )dx 1 = j bWdl  . 

J r 

x“(i 


(91) 


In  order  to  expand  dF/dx  as  given  in  (87)  in  a power 
series  as  was  done  in  the  case  in  which  £ and  b were  constants, 
it  will  be  convenient  to  write 
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(92) 


C(x)  = r3\(x) 
b(x)  = bji(x) 

where  t and  b are  constants,  having  the  dimensions  and  order  of 
magnitude  of  C(x)  and  b(x),  respectively,  X(x)  and  p(x)  are 
thus  non-dimensional  and  of  order  one.  With  these  substitutions 
(90)  and  (91)  become 

Jx 


C 

T> 


jj. (-t)d-t 


pX+^?  |'X 

^XCx’Mx1  +j  X(xf)dx’ 


r 


0 

■x 

X 


0 


"*1 


ji(-t)d  l 


f'x^r,  ,>x 

i 1 X(x’)dx‘  + X(x,)dxt 


J 0 


J 0 


(o0n 


(91') 


In  (90l)  and  (91f)  we  now  consider 

- C(*) 
e b(x) 


(92') 


as  a function  of  ^ and  ^ respectively,  from  which  we  finally 
obtain  power  series  in  e for  and  Thus  from  (90') 


e(0)  = 0 


.aicai  = 


(x) 


ad 


2 i \(x 1 )dx ' 

JO 

p'(x)  ; ^x'Mx'  - u(x)  X(x) 
v'O , 

r 1 x 2 

2[  v X(x* )dx 1 ] 
d 0 


(93) 


from  which,  considering  as  a function  of  e,  we  have 
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^2(o) 


d?;2(c) 

d e 


d%<0) 


a £ 


= 0 


fx 

2 X(x’)dxf 
^0 

li(x) 


*+  \ X(x')dx' 
!0 


/(*) 


X(x*  )dx 1 


Similarly,  from  (91')  we  have 
e(0)  = 0 


ddOi  = 
d£n 


p>(x) 
tx — 


2 J X(x»)dx' 


d2e(0) 


-uf(x)  f X(xf)dxr'  + y,(x)X(x) 

_ J 0 - 


rx  p 

2[  X(x’  )dx'  1 


0 


from  which,  considering  ^ as  a function  of  e. 


^(0) 


d$  (0) 
dc  ’ 


= 0 


2 I XCx^dx’ 

J 0 


p.(x) 


„ rx 

dtfJL<0)  ^ I X(x f )dx 1 
1 _ vO 


de£ 


,3(x) 


[ u’(x)  \(x ’ )dx 1 - 

6 


(9h-) 

- M-(x)  X(x)  1 « 

(95) 

we  have 

(96) 

^i(x)  XCx)  3 • 


Using  (92)  we  may  write  (87)  in  the  form 
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a(x)[p.(x)  ~eX(x) 1 


dF 

dx  a £si  +p,  ( x ) — e\(x 
b 


- fS2(X’  2S  + 0)  - S2tX  - f’l’  ^ + °)] 

(87') 


a(x+£0)[n(x)  +eX(x) ] r , n v „ rtvi 

[s  (x,  ■*-  -1-  0)-s  (x  + -x-  + Q)J. 

. „ > - 2 7 2u)  2 27  2u 

a(x+£~) 

+ ^i(x-f^2)-eX:x  + j;2) 

Using  (87'),  ( 9 4-)  and  (96).  we  now  expand  dF/dx  in  a power  series 
in  e,  keeping  only  the  terms  of  lowest  order  in  £ and  obtain 

x -.2 


dF  _ }+  jy,  j~ a(x)»  b(x)  i J 0 )dX)  QSgCx, 

dx  " ax  ^aCx)+b(x)  \ b(x)  ' ax' 


+ 0) 


(97) 


in  which  we  have  substituted  (92)  and  ( 92  * ) in  order  to  return 
to  the  variables  C(x)  and  b(x).  Comparing  (97)  with  (36)  and 
following  tne  methods  used  to  arrive  at  (4-2)  and  (4-3),  we Tsee 
that  the  salinity  formerly  given  by  (b-2)  is  now  expressed  by  the 
following  similar  equation: 

-(  k(x  ’ )dx 1 
OL 

s(x)  = s(L)e  (98) 

where 

k(x) 

which  is  identical  with  (4-2)  if  aCx),  b(x)  and  £(x)  are  constant. 

Extension  of- Analysis  to  Hare  General  Kinetic  Conditions. 

Finally,  we  consider  a model  which  is  identical  to  the 
first  model  (a.  b,  % constant,  complete  mixing  at  high  and  low 
tide  and  no  mixing-  at  other  times)  expect  that  the  bottom  layer 
is  assumed  to  have  a periodic  velocity 


2iB  a(x) 

2u)  a(x) 


t b(x) 
* b(x) 


h 


b(x) 


i2 


(98*) 


C(x»)dx*  | 

0 J 
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u-^(t)  = - U cos  ( cat  + 9). 


39 

(99) 


For  this  case  the  flux  formerly  given  hy  (28)  is  now 

5’*V’2w-0 

f 

F(x)  = | (b9u2s2  + b1u1s1)dt 

Jis/2u-0 


3 TT/  2 'it-  0 


h/2uh-0 


?*/2« 


f(x,t)dt  + | 


f(x,t)dt 


J 


3n/2 


(100) 


where 

f(x,t)  =■  b2u2 s2  + biii^i,  b2  = b + C sin  oJt,  b-^  = a 

and  u2  i^  given  by  (6),  and  Dy  (99)» 

Following  the  work  done  previously  we  now  have,  in 
place  of  (29), 


since 


4E  ~ (b  + C)s2(x,  JL  + 0)  + a S-,  (x,  -2L  + 0) 

dx  ^ 2u  1 2 u> 

- (b  + Os2(x,  ~ - 0)  - a sx(x,  |S  -.0) 
from  the  continuity  Eq,  (32), 


(101) 


(bisi  + bpSo)j  = (b-iSn  + b9Sp}  , 

- 13ti/2u-0  1 ‘ ^ ^'3n/2u)^0 

Since  we  have  assumed  mixing  complete  after  high  and  low  tides, 

(15)  and  (16)  may  still  be  used.  The  continuity  equations  (31) 

and  (32)  may  also  be  applied  to  this  model  and  since  u2  is  still 

given  by  (6),  Eqs.  (9),  (10),  (11),  and  (12)  may  also  be  used. 

However,  (13)  and  (lh)  must  now  be  changed  to  (see  Eqs.  (9)  and 

(ID) 
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s-|(x, 


(4n  + Dr  „ an  _ „ . (4n  - Pit  + C) 

'2*  2 w 


2u 


- 0)  = s^(x  + 


S (X,  - 0)  = s,Xx  - Yin 

1 2u 

n = 0,  1 , 2 , • » • 


where 


(4n  + l)rc 

2«J 


0) 


p(4n+l)n/2w 


Yc 


u 


u (t)dt  = ~ -S-2XLJE. 

1 U) 


(4n-l)u/2w 
f'(4n+3)'rc/2u 

u^(t)dt  = 
(4n+l)Tc/2w 


2U-  cos 


2 


w 


(102) 

(103) 


( 102  * ) 


(103*) 


Since  with  b^  = a,  a constant,  we  have  y?  = Y]>  we  denote  their 
common  value  by 


20  cos  © 

r - — - > 


(104) 


u 


the  distance  traveled  by  the  lower  layer  between  high  and  low 
tide. 

v/e  now  express  the  various  terms  in  (101)  in  terms  of 
Sp(x.  —■  + 0)-. 


From  (15)  s-,(x,  JL  + 0)  = s?(x,  ~ + 0)  and 

1 2u  d 


2w 


from  (16)  s^(x,  -^2  + o)  = s2(x,  ~ + 0)  and 

from  (11)  s2(x,  - 0)  = s2(x  - Kv  ^ + 0) 

where  x - is  given  in  (12). 

From  (103),  (15)  and  (104)  we  have 

it 


(x,  4-  " 0)  = Sp(x  - y,  — + 0), 


2cj 


<ZU) 


(105) 


From  (16) 


s (x 

SliX’  2^ 


:)  — s p ( x , 4"  + 0) 
2u> 
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and  substituting  (11),  (1C5)  and  (16)  in  (32)  we  have 

S2(x,  |s  + 0)  = [a  s?(x-Y,  JE.  40)+<b-C)s2(*-?1J  JL  +0)], 


fl  r/.  > 


From  (9)  ana  (106) 


S2(X>  Is  - 0)  = Ilfex  U SS(X+61>5S  + 0)+(b-C)s2(x,  + 0)3 


where 


61  = ^2  “ Y 


(107) 

(108) 


following  the  method  used  for  evaluating  y» 

Finally,  from  (102),  Xl6)  and  (106) 

sT(x,  o)  = — J— p [a  sP(x,  -1  + 0)+(b-C)s?(x-6  , JL  + o)  ] 

1 2u  a+b-C  *=  2u  2 27  2w 

(109) 


where,  from  (7’) 


* _ (x  + Y)(b  - C) 
x " 2 b^Tc 


(110) 


Substituting  (11),  (15),  (16),  (105),  (106),  (107),  and  (109) 
in  (101)  we  have 


1 = frtTjf  ts2<x>  ■£,  + °)  - s2<x  + 61>  h + 0)] 


a+b-i 


(111) 


[s2(x,  iL  + 0)  - s2(x  - 62,  + 0)]  . 


2’ 


If  we  assume  y to  be  of  order  C/b  and  that  e = C/b  « 1 we  may 


write 


[s0(x,  ~ + 0)  - Sp ( x + 6-,,  ~ + 0)1 
dx  a+b(l-e)  *-  2w  d ~ 2u  (111* 

+ ^nri7  [s2U’  -fc  + 0)  ■ S2(x  ~ &2»  *■  °>  3 


where  from  (10),  (12),  (108),  and  (110) 


2xe 

6i = rrr  - r 


f 11  O'! 


tff.  /O 
lioriA  - £.\j£.f  d. 


^ * 1 " ■'  ■' 
2 1 4 £ 


0.13) 


Using  (1111),  (112),  (llj),  and  (IC'-f)  we  expand  the  right  hand 
side  of  (ill’)  in  a power  series  in  e and  keep  only  the  terras 
of  lowest  order  in  e = C/b,  The  result  is 


n ^ „ 3So(x,  + 0) 

df  _ = ab  9 r / _v  \ 2 1 

dx  a+b  3x  b ax 


(114) 


This  result  is  just  what  one  might  have  expected  since,  to  first 
order  in  e,  the  distance  traveled  by  the  upper  layer  in  going 

Oy  y 

from  low  to  high  tide  is,  as  was  previously  shown,  £ = — ~ , 

whereas  the  distance  now  traveled  by  the  lower  layer  in  this 
■ - - 1 -i  • . 2 J cos  cp  , , ii j . .u < . u 


same  time  interval  is  y = 


; and  hence  the  term  which 


corresponds  to  tiie  diffusion  coefficient  is  again  proportional 
to  the  square  of  the  relative  displacements  of  the  top  and  bottom 
layers.  Following  the  work  previously  done,  we  could  also  solve 
(114)  and  obtain  the  salinity.  The  result  is 


s(x)  = sXDe 


L-Xn ' X-X, 


where 


ix0s 


- Y * 


(115) 


mati.cn  in 


Non- steady  Case  s 


Presence 


of  a Sc  lute  Source  Distribution. 

Let  us  return  to  equation  (4l).  The  left  hand  side 
expresses  the  flux  of  salt  upstream  pact  the  cross-section  at  x 
during  a tidal  cycle  as  a consequence  of  the  'tides.  The  right 
hand  side  expresses  the  flux  downstream  past  the  same  section 
and  during  the  same  interval  due  to  the  river.  Equation  (4l) 
expresses  the  fact  that  these  two  contributions  to  the  flux:  must 
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cancel  in  an  estuary  in  which  conditions  do  not  change  from  cycle 
to  cycle  and  in  which  no  salt  (or  other  solute)  is  introduced 
into,  or  removed  from,  the  estuary  by  external  agents. 

We  now  allow  for  variation  of  the  solute  concentration 
from  cycle  to  cycle  and  also  for  introduction  of  solute  into 
the  estuary. 


s(x+Ax,  t) 
x+Ax 

Consider  an  element  of  estuary  of  length  Ax  and  of  mean 
c-rooS- sectional  area  during  a tidal  cycle  a+b.  The  flux  out  of 
the  element  at  x during  the  time  At  due  to  the  tides  is: 


^ ab_  “■  x2 
"~^2  a+b  X 


-2LS  At 


The  flux  into  the  element  at  x + Ax  during  time  At 
due  to  the  tides  is  (neglecting  terms  which  will  vanish  in  the 
limit  Ax,  At— >0) 


At  is 


=5  5b  5,  tx2  ||  ♦ £ (X2  ff  )Ax]At  - g At. 


bi  a+b  2*  lx  S +-S  u ax-  - 
The  net  flux  into  the  element  due  to  the  tides  during 


4C_ 

b“ 


■Sfe ..  — hf  _2-  (x^  •^)Ax/vt 
a+b  2 it  dx  ax;AxaT 


t,AXAt- 


(116) 


Let  us  assume  that  the  derivation  of  expression  (11 6)  remains 
valid  if  we  permit  changes  in  solute  concentration  from  cycle 
to  cycle  as  well  as  solute  introduction  by  external  agents 
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(provided  our  fundamental  mixing  model  is  net  altered).  The  tine 
deoendence  of  the  salinity  will  most  likely  be  due  to  variation 
in  the  river  discharge.  The  external  agent  may  be  the  discharge 
of  pollution  by  an  industry.  The  mixing  model  may  then  be  ex- 
pected to  be  essentially  unchanged  provided  the  time  during  which 
the  river  discharge  varies  significantly  is  large  compared  to  a 
tidal  cycle  and  provided  the  amount  of  solute  which  is  introduced 
during  a tidal  cycle  is  so  small  that  it  does  not  materially 
affect  the  concentration  during  this  cycle. 

The  flux  into  the  element  in  time  At  at  x due  to  the 


river  is 


R( t)sAt. 


is 


The  flux  out  of  the  element  at  x + Ax  due  to  the  river 


R(  s + M Ax)At. 
ax 


The  net  flux  out  of  the  element  due  to  the  river _ is 


R(t)  AxAt. 

ox 


If  q(x,t)  d enntes  the  introduction  of  solute-  into  the 
estuary  p^r  unit  length  of  ostuary  and  per  unit  time  (dimensions 
ML  ^T”^ ) than  the  total  increase  of  solute  in  the  clement  during 


At  is 


R (x2  ~*)AxAt  ~ R AxAt  + q(x,t)AxAt. 


dx 


(118) 


This  must  be  equal  to 

[s(t  + At)  - s(t)  ] (a  + b)Ax 


(a  + b)AxAt« 


(119' 


»T«  T*C.  fes 
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Hence  we  obtain  the  following  differential  equation  for  the 
concentration  of  solute 


( aj-b) 


at 


+ R(t) 


££  = 
ax 


r> 


ax 


(x2  |5) 
ax 


+ q (x , t ) . 


v J.xy; 


Conclusion* 

In  this  report  we  have  derived  an  equation  governing 
the  distribution  of  a solute  in  an  estuary  in  which  the  river 
flow  and  the  tides  are  the  predominant  factors  and  in  which  the 
dynamic  and  mixing  conditions  can  be  described  by  a very  simple 
physical  model.  This  model  was  subsequently  generalized  in 
various- 'ways  and  it  was  found  that  the  basic  form  of  the  diffu- 
sion equation  was  unchanged*  It  should  be  noted  that  the  various 
generalisations  which  were  discussed  independently  can  all  be 
combined  into  one  model.  In  the  last  section  an  equation  was 
derived_which  allows  for  time  dependence  of  the  various  quanti- 
ties involved  and  for  introduction  of  solute  Into  the  estuary  by 
an  external  agent* 

In  order  to  ascertain  if  the  theory  is  capacie  or  des- 
cribing the  conditions  in  a real  estuary  it  will  be  necessary  to 
compare  the  predictions  of  this  theory  with  the  observational 
data  from  a number  of  estuaries  whose  over-all  structure  permits 
application  of  this  theory.  In  particular  this  will  imply  that 
the  estuaries  are  essentially  vertically  homogeneous,  that  their 
properties  arc  determined  mainly  by  the  river  and  tides  and  that 
it  is  feasible  to  divide  the  cross-section  at  any  position  along 
thu  estuary  into  two  regions  such  that  the  tidal  velocity  is 


\T 

i 


o nr- 
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essentially  uniform  within  each  region  but  different  for  the  two 
regions. 


However,  even  if  it  is  found  that  this  theory  is  net 
directly  applicable  to  real  estuaries  it  may  nevertheless  be  a 
helpful  guide  in  experimental  studies.  Moreover,  it  is  especi- 
ally to  be  hoped  that  it  will  be  a useful  bfcsis  for  the  formula- 
tion of  more  elaborate  theories  that  will  describe  real  estuaries* 
A few  interesting  solutions  of  the  diffusion  equation 
derived  in  the  last  section  will  be  given  in  a future  report. 
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